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Cubic  Spline  Interpolation  on 
Nested  Polygon  Triangulations 


Oleg  Davydov,  Gunther  Niirnberger,  and  Frank  Zeilfelder 


Abstract.  We  develop  an  algorithm  for  constructing  Lagrange  and  Her- 
mite  interpolation  sets  for  spaces  of  cubic  C^-splines  on  general  classes  of 
triangulations  built  up  of  nested  polygons  whose  vertices  are  connected 
by  line  segments.  Additional  assumptions  on  the  triangulation  are  signifi- 
cantly reduced  compared  to  the  special  class  given  in  [41.  Simultaneously, 
we  have  to  determine  the  dimension  of  these  spaces,  which  is  not  known 
in  general.  We  also  discuss  the  numerical  aspects  of  the  method. 


§1.  Introduction 

In  contrast  to  univariate  splines,  it  is  a non-trivial  problem  to  construct  even 
one  single  set  of  interpolation  points  for  bivariate  spline  spaces.  Such  interpo- 
lation sets  for  Sq( A),  the  space  of  splines  of  degree  q and  smoothness  r,  were 
constructed  for  crosscut-partitions  A (see  the  survey  [9]  and  the  references 
therein).  For  general  triangulations  A,  interpolation  sets  were  constructed  for 
51(A),  q > 4 in  [3]. 

The  case  g = 3 is  much  more  complicated  given  that  not  even  the  dimen- 
sion of  5^(A)  is  known  for  arbitrary  triangulations  A.  It  is  an  open  question 
whether  the  dimension  of  5]  (A)  is  equal  to  Schumaker’s  lower  bound  [12],  The 
aim  of  this  paper  is  to  investigate  interpolation  by  5^  (A)  for  general  classes  of 
triangulations  A consisting  of  nested  polygons  whose  vertices  Eire  connected 
by  line  segments.  Following  a general  principle  of  locally  choosing  interpo- 
lation points  for  53(A)  by  passing  from  triangle  to  triangle,  we  describe  an 
inductive  method  for  constructing  point  sets  that  admit  unique  Lagrange  (re- 
spectively Hermite)  interpolation  by  5|( A)  under  certain  assumptions  on  A. 
Moreover,  we  prove  that  the  dimension  of  these  spaces  is  equal  to  Schumaker’s 
lower  bound. 

In  this  way  we  obtain  a class  of  triangulations  A which  is  significantly 
larger  than  the  special  class  described  in  [4].  Moreover,  the  methods  of  proof 
in  this  paper  are  different  from  those  in  [4].  It  is  important  to  note  that 
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triangulations  of  this  type  can  be  constructed  starting  from  any  given  points 
in  the  plane,  see  [11], 

The  numerical  examples  (with  up  to  100,000  interpolation  points)  show 
that  in  order  to  obtain  good  approximations,  it  is  desirable  to  subdivide  some 
of  the  triangles.  Our  method  of  constructing  interpolation  points  also  works 
for  these  modified  triangulations. 

We  note  that  our  interpolation  method  can  be  used  for  the  construction 
of  smooth  surfaces  without  involving  any  derivative  data.  For  scattered  data 
fitting,  the  needed  Lagrange  data  are  approximately  computed  by  local  meth- 
ods. In  contrast  to  the  finite  element  methods  for  cubic  splines,  we  do  not 
need  to  subdivide  all  triangles  by  a Clough- Tocher  split  or  use  derivatives. 


§2.  Preliminaries 

Let  A be  a regular  triangulation  of  a simply  connected  polygonal  domain 
n in  H2.  We  denote  by  53(A)  = {s  £ C'1(fl)  : s|t  € II3,  T £ A}  the 
space  of  bivariate  splines  of  degree  3 and  smoothness  1 (with  respect  to  A). 
Here  n3  = span  {xvy>1  : v,  p > 0,  v + p < 3}  denotes  the  space  of  bivariate 
polynomials  of  total  degree  3. 

We  investigate  the  following  interpolation  problem.  Construct  a set 
{zi, . . . , 2 at}  in  SI,  where  N — dim  S3  = (A),  such  that  for  each  function  / 6 
C(Sl),  a unique  spline  s € S3  (A)  exists  such  that  s(z,)  = f(z ,),  i = 1, . ..,N. 
Such  a set  (zj, . . . , zjv}  is  called  a Lagrange  interpolation  set  for  S-j(A).  If  also 
partial  derivatives  of  / are  involved,  then  we  speak  of  a Hermite  interpolation 
set  for  S£(A). 

In  contrast  to  [4],  we  will  use  Bernstein- Bezier  techniques  [2,5].  Given  a 
spline  s € S^A),  we  consider  the  following  representation  of  the  polynomial 
pieces  p = s|r  € n3  on  the  triangle  T € A with  vertices  vi,v2,  V3, 


p{x,y) 


53  (x,y)mx,y)mx,y),  (x,y)€T, 

V+p+<7=3 


(1) 


where  $1  € l = 1,2,3,  is  uniquely  defined  by  <t>i(vk)  = 6k, 1,  k = 1,2,3. 
This  representation  of  p is  called  the  Bernstein- Bezier  representation  of  p, 
the  real  numbers  a^l<T  are  called  the  Bernstein-Bezier  coefficients  of  p,  and 
$i(x,y),  l = 1,2,3,  are  the  barycentric  coordinates  (w.r.t.  T)  of  ( x,y ) € T. 

Definition  1.  A set  A C {(v,  p,  a,  T)  : v + p + <7  = 3,  T e A}  is  called  an 
admissible  set  for  S^(A)  if  for  every  choice  of  coefficients  , (is,  p,o,T)  £ A, 

a unique  spline  s € S3  (A)  exists  with  these  coefficients  in  the  above  Bernstein- 
Bezier  representation. 


The  above  Bernstein-Bezier  form  can  be  used  to  express  smoothness  con- 
ditions of  polynomial  pieces  on  adjacent  triangles  T\ , T2  with  vertices  rq , i>2,  V3, 
respectively  wi, U2,  U4  (cf.  [2,5]). 


Theorem  2.  Let  s be  a piecewise  cubic  polynomial  function  defined  on  Ti  U 


T2.  Then  s 6 S^({TUT2})  iff  a™0  = a'T‘|0,  u + p = 3, 


aSL,o$i(u4)  + v + p.  = 2. 


[A] 


and  aS!i 
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For  later  use,  we  also  mention  here  the  following  relations  between  the 
Bernstein-Bezier  coefficients  of  a cubic  polynomial  p in  the  representation  (1) 
and  its  partial  derivatives  at  v\  in  direction  of  a unit  vector  parallel  to  the 
edge  e = [vi,v2\,  denoted  by 


a' 


[T] 

3,0,0 


:pM, 


a: 


[T] 

2,1,0 


V2W2 


al,2,0 

dpjvi) 

de 


3 de 
\t] 


[T] 


■Ml 


(2) 


3(a2,i,o  — a3,o,o)  d2p(v\)  _ 6(a^2  0 p + 03^0) 


,P1 


bl 


■ V2\\2 


de 2 


IK  - V2 


§3.  Main  Results 

In  this  section,  we  state  our  main  results  on  £3 (A),  where  A consists  of 
nested  polygons  whose  vertices  are  connected  by  line  segments.  We  first  define 
this  class  of  triangulations.  Then,  we  determine  the  dimension  and  construct 
interpolation  sets  for  the  corresponding  spline  space.  Moreover,  we  show  that 
this  dimension  is  equal  to  Schumaker’s  lower  bound  [12].  Finally,  we  discuss 
a property  of  A which  is  essential  for  the  local  construction  of  interpolation 
points. 

First,  we  describe  triangulations  of  nested  polygons  and  decompose  the 
domain  into  finitely  many  subsets  needed  in  our  construction  of  interpolation 
points. 

Triangulations  of  nested  polygons.  We  consider  the  following  general  type 
of  triangulation  A.  Let  Po,  Pi, . . . , Pk  be  a sequence  of  closed  simple  polygonal 
lines,  and  let  be  the  closed  (not  necessarily  convex)  bounded  polygon  with 
boundary  P^.  Suppose  that  the  polygons  are  nested,  i.e.,  C 

p = 0, . . . , k.  The  vertices  of  A are  the  vertices  of  PM,  p = 0, . . . , k,  and 
one  vertex  inside  Po ■ The  edges  of  A are  the  edges  of  P p = 0, . . . , k,  and 
additional  line  segments  connecting  the  vertices  of  P ^ with  the  vertices  of 
Pfi+\,  p = 0, . . . , k — 1.  The  resulting  triangulation  A of  O :=  does  not 
have  vertices  in  the  interior  of  \ p = 0, . . . , k — 1,  and  does  not  have 
edges  connecting  two  vertices  of  P^  other  than  the  edges  of  PM,  see  Figure  1. 
Decomposition  of  the  domain.  We  decompose  the  domain  into  finitely 
many  sets  Vo  C Vi  C • ■ • C Vm  = Q,  where  each  set  Vi,  is  the  union  of  closed 
triangles  of  A,  i = 0, . . . ,m.  Let  Vo  be  an  arbitrary  closed  triangle  of  A in 
flo-  We  define  the  sets  Vi  C • • • C Vm  by  induction.  Assuming  V)_  1 is  defined, 
we  choose  a vertex  vt  of  A such  that  there  exists  at  least  one  triangle  of  A 
with  vertex  w,  and  a common  edge  with  Vi- 1.  Let  T^i, . . . ni  > 1, 

be  all  such  triangles.  We  set  Vi  = Vj_i  U Thi  U • • • U TitTl% , and  denote  by 
A i = {T  g A : T C K}  the  subtriangulation  which  corresponds  to  the  set  V,. 

The  vertices  Vi,  i = 1, ...  ,m,  are  chosen  as  follows.  After  choosing  Vo  to 
be  an  arbitrary  closed  triangle  of  A in  fio>  we  pass  through  the  vertices  of  Po 
in  clockwise  order  by  applying  the  above  rule.  (It  is  clear  that  the  choice  of 
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Fig.  1.  Triangulation  of  nested  polygons. 

these  vertices  is  unique  after  fixing  the  first  vertex.)  Now,  we  assume  that  we 
have  passed  through  the  vertices  of  PM_],  We  fix  a vertex  w M of  P^  that  is 
connected  with  at  least  two  vertices  of  PM+i-  Then  w.r.t.  clockwise  order,  we 
choose  the  first  vertex  of  PM  greater  than  which  is  connected  with  at  least 
two  vertices  of  P,,  _ i . Then  we  pass  through  the  vertices  of  P(J  in  clockwise 
order  until  w~,  and  pass  through  the  vertices  of  P,,  in  anticlockwise  order 
until  by  applying  the  above  rule.  (Here  u>£  denotes  the  vertex  next  to 
in  clockwise  order  and  w~  denotes  the  vertex  next  to  in  anticlockwise 
order.)  Finally,  we  choose  the  vertex  wfl.  (It  is  clear  that  after  fixing  w^,  the 
choice  of  the  vertices  on  PM  is  unique.) 

The  construction  of  an  admissible  set  for  £3  (A)  and  the  choice  of  inter- 
polation points  depend  on  the  following  properties  of  the  triangulation  A. 

Definition  3.  (1)  An  interior  edge  e with  vertex  v of  the  triangulation  A is 
called  degenerate  at  v if  the  edges  with  vertex  v adjacent  to  e lie  on  a line. 
(2)  An  interior  vertex  v of  A is  called  singular  if  v is  a vertex  of  exactly  four 
edges  and  these  edges  lie  on  two  lines.  (3)  An  interior  vertex  v of  A on  the 
boundary  of  a given  subtriangulation  A'  of  A is  called  semi-singular  of  type 
1 w.r.t.  A'  if  exactly  one  edge  with  endpoint  v is  not  contained  in  A'  and 
this  edge  is  degenerate  at  v.  (4)  An  interior  vertex  v of  A on  the  boundary 
of  a given  subtriangulation  A'  of  A is  called  semi-singular  of  type  2 w.r.t.  A' 
if  exactly  two  edges  with  endpoint  v are  not  contained  in  A'  and  these  edges 
are  degenerate  at  v.  (5)  A vertex  v of  A is  called  semi-singular  w.r.t.  A'  if  v 
satisfies  (3)  or  (4). 

In  the  following,  we  construct  an  admissible  set  and  interpolation  sets  for 
1S3  ( A),  where  A is  a nested-polygon  triangulation. 

Construction  of  an  admissible  set.  First,  we  choose  Aq  = {(v,  p,  a,  Vo)  : 
v+p+cr  = 3}  and  then,  proceeding  by  induction,  we  successively  add  admissi- 
ble points  on  V,  \ V)_i,  i = 1 ,...  ,m.  Assuming  that  an  admissible  set  A,-i  on 
Vi- 1 has  been  constructed,  we  choose  admissible  points  on  V,  \ V, _ 1 as  follows. 
By  the  above  decomposition  of  ft,  V*  \ is  the  union  of  consecutive  trian- 
gles Tip, . . . ,Ti>ni  with  vertex  v,  and  common  edges  with  V) _ 1 . We  denote 
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the  consecutive  endpoints  of  these  edges  by  ■ ■ • , , and  the  piece- 

wise  polynomials  in  the  representation  (1)  on  T,j  by  pij  G II3,  j = 1, . . . nt, 
where  the  vertices  of  Tij  are  ordered  as  follows:  Vi,Vij,Vij+±.  Furthermore, 
we  denote  by  eij  the  edges  [vij,Vi],  j — 0, . . . ,n*. 

We  need  the  following  properties  of  the  subtriangulation  A*  = {T  G A : 
T C Vi}  at  the  vertices  Vit 0, . . . , 

(a)  eij  is  non- degenerate  at  Vij, 

(b)  Vij  is  semi-singular  w.r.t.  A;.  (This  latter  property  is  only  relevant  if 
Vij  lies  on  the  boundary  of  A*,  i.e.,  for  j G {0,71;}.) 

For  j G {1, . . . , n<  — 1},  we  set  Cij  = 1 if  (a)  holds,  and  c,j  = 0 otherwise.  For 
j G {0,7i;},  we  set  c;  j = 1 if  both  (a)  and  (b)  hold,  and  c;j  = 0 otherwise. 
Moreover,  we  set  c;  = Yl’jL 0 c«,j)  and  assume  c;  < 3,  i = 1, . . . , 71. 

Now,  we  construct  the  following  admissible  points  on  V;  \ V;_i.  If  c;  = 3, 
then  no  point  is  chosen.  If  c;  = 2,  then  we  choose  (3, 0, 0, 7}i).  If  c;  = 1,  then 
we  choose  (3,0, 0, T^i)  and  (2, 0, 1, Tij),  where  e^j  is  an  edge  with  c;j  = 0. 
If  c;  = 0,  then  we  choose  (3,0,0,T;ii),(2,0, 1,7;^)  and  (2,  l,0,T;ti).  The 
admissible  set  Ai  on  Vi  is  obtained  by  adding  these  points  to  Ai-\. 

Construction  of  interpolation  sets.  We  choose  interpolation  points  in  Vo 
and  then  in  V;  \ F;_i,  i = 1, . . . ,7n,  successively.  In  the  first  step,  we  choose  10 
different  points  in  Vo  (respectively  10  Hermite  interpolation  conditions)  which 
admit  unique  Lagrange  interpolation  (respectively  Hermite  interpolation)  by 
the  space  n3.  For  example,  for  Lagrange  interpolation,  we  may  choose  four 
parallel  line  segments  in  Vo  and  v different  points  on  each  lv,  v = 1, 2, 3, 4. 
Assuming  that  the  interpolation  points  in  V;_i  have  already  been  chosen,  we 
proceed  to  V,  \ V;_i  as  follows. 

For  Lagrange  interpolation,  we  choose  the  following  points  in  V;  \ V;_i. 
If  c;  = 3,  then  no  point  is  chosen.  If  c;  = 2,  then  we  choose  u;.  If  c;  = 1,  then 
we  choose  u;  and  one  further  point  on  some  edge  eij  with  Cij  = 0.  If  c;  = 0, 
then  we  choose  u;  and  two  further  points  on  two  different  edges. 

For  Hermite  interpolation,  we  require  the  following  interpolation  condi- 
tions for  s G 53(A)  at  the  vertex  i>;.  If  c;  = 3,  then  no  interpolation  condition 
is  required  at  v;.  If  c;  = 2,  then  we  require  s(u;)  = /(u;).  If  c;  = 1,  then 
we  require  s(w;)  = /(v;)  and  gf£-(u;)  = g^-(v;),  where  eij  is  some  edge 

with  Cij  = 0.  If  c;  = 0,  then  we  require  s{vi)  = f(vi),  §§(vi)  = g^{vi)  and 

&(«<)  = &(«!). 

By  the  above  construction,  we  obtain  a set  of  points  for  Lagrange  inter- 
polation respectively  a set  of  Hermite  interpolation  conditions. 

Theorem  4.  Let  A be  a triangulation  of  nested  polygons.  If  for  all  i G 
{1, . . . ,771},  c;  < 3 and  no  vertex  u;  is  simultanously  semi- singular  (of  type 
2)  w.r.t  A;  and  non-singular,  then  a unique  spline  in  5} (A)  exists  which 
satisfies  the  above  Lagrange  ( respectively  Hermite)  interpolation  conditions. 
In  particular,  the  total  number  of  interpolation  conditions  is  equal  to  the 
dimension  of  S}(A). 
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Proof:  First,  we  prove  that  the  set  constructed  above  is  an  admissible  set 
for  53(A).  To  this  end,  we  show  by  induction  that  A,  is  an  admissible  set  for 
5!(A)|Ai  - {s|Ai  : s £ 53(A)}.  This  is  clear  for  i = 0.  Now,  we  assume  that 
Ai-i  is  an  admissible  set  for  53(A)|Aj_,,  where  i £ {1, . . . ,m},  and  consider 
Vi-  For  simplicity,  we  omit  here  the  index  i for  vt,  vt  j,  ej.y , Pjj , 7}.-  and  7ij. 
It  follows  from  the  induction  hypothesis  and  Theorem  2 that  the  coefficients 
0^3 a = 0, . . . ,3  - v,  v = 0, 1,  of  p,  £ II3,  j = 1, . . . , n,  on  Tjt  are 
uniquely  determined.  Moreover,  if  Cij  = 1 for  some  j £ {1, . . . ,n  — 1},  then 

[T  1 

it  follows  from  Theorem  2 that  the  coefficient  a'2  is  uniquely  determined. 

In  the  following,  we  show  that  if  c,  o = 1,  then  the  coefficient  o^W,  is 
uniquely  determined.  Let  us  consider  the  case  where  vq  is  semi-singular  of 
type  2 w.r.t.  A,.  (The  case  that  vq  is  semi-singular  of  type  1 w.r.t.  A;  is 
analogous.)  We  denote  by  T;  £ A,  l = 1,...,3,  the  triangles  with  vertex  Vq 
not  contained  in  A,  in  anticlockwise  order,  and  by  e;  the  common  edge  of  7} 
and  T;+i,  l = 1,2.  Since  T3  has  a common  edge  with  A;_i,  it  follows  from 

Theorem  2 that  the  coefficient  of  p3  £ IT3  on  T3  is  uniquely  determined. 
Moreover,  since  £2  and  ei  are  degenerate  at  v0,  the  coefficients  of  pi  £ II3 
on  7),  Z = 1,2,  are  uniquely  determined.  Since  eo  is  non-degenerate  at  Vo,  it 
follows  from  Theorem  2 that  the  coefficient  a^i.o 's  uniquely  determined.  We 
note  that  since  A is  a nested-polygon  triangulation,  at  least  two  triangles  with 
vertex  vo  not  contained  in  A,  exist.  Therefore,  if  c^o  = 0,  then  the  coefficient 
is  not  yet  determined. 

Analogously  as  above,  it  can  be  shown  that  the  coefficient  a2Q  \ is  uni- 
quely determined  if  Cj,n  = 1.  Otherwise,  this  coefficient  is  not  yet  determined. 

Now,  we  consider  the  vertex  v.  The  arguments  below  will  show  that  we 
may  assume  that  v is  an  interior  point  of  A.  We  denote  by  Tn+i  £ A,  l = 
1, . . . ,r,  r > 3,  the  triangles  with  vertex  v not  contained  in  A*  in  anticlock- 
wise order.  Moreover,  let  the  piecewise  polynomials  pn+i  £ II3,  l = 1, . . . , r, 
on  Tn+i  in  the  representation  (1)  be  given  such  that  the  first  barycentric  co- 
ordinate always  corresponds  to  v.  The  above  arguments  show  that  exactly 
Cj  < 3 coefficients  of  the  set  C\  — a : cr  = 0, . . . , 3 — u,  v — 

2,3,  1 = l,...,n  + r}  are  uniquely  determined.  On  the  other  hand,  we  con- 
struct 3 — Cj  additional  admissible  points  from  C\  on  V;  \ V,  _ 1 . Now,  it  follows 
from  the  Cll-property  at  v and  Theorem  2 that  all  coefficients  from  Ci  are 
uniquely  determined.  By  our  method  of  passing  through  the  vertices  of  A, 
v is  not  semi-singular  of  type  1 w.r.t.  Aj.  In  particular,  if  v = wfl  for  some 
//  £ {0,  Moreover,  by  assumption  v can  be  semi-singular  of  type  2 

w.r.t.  Aj  only  if  v is  singular.  In  this  case,  we  have  r — 3,  and  it  follows 
from  Theorem  3.3  in  [13]  that  the  coefficient  a^j  is  uniquely  determined. 
Otherwise,  if  r > 4,  then  for  some  l £ {1,  ...,r  — 1}  one  common  edge  of 
Tn+i  and  Tn+(+1  is  non-degenerate  at  v,  and  we  can  also  proceed  with  our 
arguments. 

Since  all  relevant  differentiability  conditions  at  the  edges  with  endpoint 
v,  respectively  Vj,  were  involved,  the  above  shows  that  Ai  is  an  admissible  set 
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for  S^AJIa*-  Thus,  the  set  Am  is  an  admissible  set  for  53(A). 

Therefore,  the  cardinality  of  Am  is  equal  to  the  dimension  of  53(A).  By 
construction,  it  is  evident  that  the  number  of  Lagrange  interpolation  points, 
respectively  the  number  of  Hermite  interpolation  conditions  coincides  with 
this  cardinality. 

By  an  inductive  argument,  it  follows  from  (2)  that  the  Hermite  inter- 
polation conditions  at  v determine  the  Bernstein-Bezier  coefficients  of  the 
admissible  points  chosen  on  V*  \ VAi-  Analogously,  the  Lagrange  interpola- 
tion conditions  uniquely  determine  the  interpolating  spline  on  the  edges  of 
Vi  \ Vi- 1.  Therefore,  the  interpolating  spline  is  uniquely  determined  on  all  of 
Vi  \ Vi-!.  This  completes  the  proof  of  Theorem  4.  □ 

For  arbitrary  triangulations,  Schumaker  [12]  gave  the  following  lower 
bound  L( A)  for  the  dimension  of  53(A), 

L(A)  — 3Vb(A)  + 2V/(A)  + <t(A)  + 1.  (3) 

Here,  VB( A)  is  the  number  of  boundary  vertices  of  A,  Vj( A)  is  the  number 
of  interior  vertices  of  A and  o(  A)  is  the  number  of  singular  vertices  of  A.  For 
bounds  on  the  dimension  of  bivariate  spline  spaces  see  also  Manni  [6]. 

Theorem  5.  If  a triangulation  A of  nested  polygons  satisfies  the  hypotheses 
of  Theorem  4,  then  the  dimension  of  5|( A)  is  equal  to  L{ A). 

Proof:  We  have  to  show  that  the  cardinality  of  Am  is  equal  to  L( A).  We 
prove  this  by  induction.  We  set  5(A0)  = 0 and  for  i G {1,  • • . m},  we  denote 
by  5(A,)  the  set  of  boundary  vertices  w of  A,  such  that  w = vito  and  c;?o  = 1 
(respectively  w = v;,n,  and  Q,n,  = 1)  for  some  l € {1,  Moreover,  let 

dj  be  the  cardinality  of  5(A,)  and  a;  be  the  cardinality  of  At.  We  will  show 
that 

L(Aj)  = Oj  + dj,  i = 0,  ...,m.  (4) 

This  is  evident  for  i = 0.  We  assume  that  L(A,_i)  = Oj_i  +d;_i  for  some  i G 
{1, . . . ,m}  and  consider  Vi.  We  have  VB(Ai)  = Ve(Ai_i)  - n,  + 2,  V/(Aj)  = 
V/(Aj_i)  + ni  — 1,  cr(Aj)  = cr(A;_i)  + 7 i,  where  7 j is  the  number  of  singular 
vertices  from  the  set  {vitj  : j = 1, . . . ,m  - 1}.  Since  a{  = aj_i  + 3 - c*, 
it  follows  from  the  induction  hypothesis  and  some  elementary  computations 
that 

L(Aj)  = a*  + d,_i  + Ci  + 7i  - n{  + 1. 

By  our  method  of  passing  through  the  vertices  of  A,  it  is  evident  that  if 
Vi,o  = 1 € 5(A,),  then  ^ 5(Aj_i).  In  the  following,  we  show  that  if 

Vi,m  € 5(Aj),  then  vi>ni  ^ 5(Aj_i).  First,  let  us  assume  that  vitTli  = vito 
for  some  Z € {1, — 1}.  If  ViiUi  is  semi-singular  of  type  2 w.r.t.  Aj,  then 
at  least  three  edges  of  A not  contained  in  A;  are  attached  to  170-  Hence, 
cifl  = 0.  If  VitTli  is  semi-singular  of  type  1 w.r.t.  A,,  then  the  edge  is 
non-degenerate  at  ViiTli,  since  Ci>ni  = 1.  Therefore,  i^o  is  not  semi-singular  of 
type  2 w.r.t.  Aj.  Again,  cj, 0 = 0 holds.  The  remaining  case  VitUi  = 17-1, ni_i, 
where  n,  = 1,  follows  by  the  same  arguments. 
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Now,  we  show  for  j (E  {1, . . . ,n»  — 1}  that  every  non-singular  vertex  vhJ 
such  that  is  degenerate  at  v,j  lies  in  <S(A,_i).  First,  we  consider  the  case 
j = 1.  Set  Vix  — Viti  and  let  eo  be  the  edge  that  connects  vtio  and  Viy.  We 
have  to  consider  two  cases. 

Case  1.  (The  vertices  and  u,_2  are  connected  by  an  edge  e.)  If  eo  is 
non-degenerate  at  V{t i then  = 1.  (In  this  case  Viti  is  semi-singular  of 

type  1 w.r.t.  Aj_i.)  Otherwise,  since  is  non-singular,  e is  non-degenerate 
at  Vi  i-  Thus,  Cj_ 2,ni_2  = 1-  (In  this  case  wtii  is  semi-singular  of  type  2 w.r.t. 
Aj_2.)  We  note  that  Vit \ is  not  semi-singular  w.r.t.  A,1+i,  since  at  least  three 
edges  of  A not  contained  in  Aj1+i  are  attached  to  v*, j. 

Case  2.  (The  vertices  and  «j_2  are  not  connected  by  an  edge.)  If  eo  is 
non-degenerate  at  tyi  then  we  also  have  = 1.  (In  this  case  i^i  is 

semi-singular  of  type  1 w.r.t.  Aj_j.)  We  note  that  v1:i  is  not  semi-singular 
w.r.t.  Aj1+i,  since  eo  is  non-degenerate  at  Vit j.  Otherwise,  let  e be  the  edge 
that  connects  Vi} i with  v^+i.  Since  is  non-singular,  e is  non-degenerate 
at  viti.  Thus,  Cjj+1,0  = 1.  (In  this  case  viti  is  semi-singular  of  type  2 w.r.t. 
Aj1+i.)  We  note  that  in  this  case  is  semi-singular  of  type  1 w.r.t.  A;_i, 
but  = 0. 

Now,  we  consider  the  remaining  case  j g {2, . . . ,m  — 1}.  Set  Vij  = vtj 
and  let  e be  the  edge  that  connects  Vij  with  v^+i.  Since  Vij  is  non-singular, 
it  follows  that  Vij  is  not  semi-singular  of  type  2 w.r.t.  to  A^ . Therefore,  e is 
non-degenerate  at  Vij.  Hence,  ctj+ii0  = 1.  (In  this  case  v,tj  is  semi-singular 
of  type  1 w.r.t.  A^+i.)  We  note  that  in  the  case  j € {2,...,n;  — 1},  by 
our  method  of  passing  through  the  vertices  of  A,  the  value  c^.- is  not 
influenced  by  the  geometrical  properties  of  A at  Vij. 

The  above  proof  now  implies  d;  = dj_i  + c,  -f  y,  — n,  -f  1,  and  therefore, 
(4)  holds.  Since  am  = 0,  we  get  L( A)  = am.  This  proves  the  theorem.  □ 

In  Theorem  4 we  assume  that  for  all  i g {1,  ...,m},  no  vertex  v,  is 
simultaneously  semi-singular  (of  type  2)  w.r.t.  A,  and  non-singular.  In  the 
following,  we  show  that  this  assumption  is  essential  for  the  local  construction 
of  interpolation  points. 

Example  6.  Let  v = V{  = (0,0),  V5  = v0  = vi:o  = (y,0),  7 < 0,  v\  = Viti  — 
(r,  mr),  r < 0,  m > 0,  V2  — = (0,5),  5 < 0,  and  set  V3  = (a,0),  a > 

0,  Vi  = (0,  /?),  fd  > 0.  Let  v be  connected  with  u3  and  Vi  and  «j_  1 be  connected 
with  vi,  l = 1, ...  ,5.  Then  v is  simultaneously  semi-singular  (of  type  2)  w.r.t. 
A,  and  non-singular.  Furthermore,  we  denote  by  T;  the  triangle  with  vertices 
v,  vi-i,vt  and  by  pi  g n3  the  polynomial  pieces  on  7),  / = 1, ...  5,  in  the  rep- 
resentation (1).  We  consider  the  set  C2  = {aj^3_„_tr|ix,  <r  = 0, . . . , 3 — u,  v = 
1, . . . , 3,  l = 1, . . . , 5}.  For  C^-splines,  it  follows  from  Theorem  3.3  in  [13]  that 
each  subset  of  C2  that  uniquely  determines  all  coefficients  of  C2  has  cardinality 
8 and  contains  the  coefficients  ofj^,  l = 3,4.  By  the  proof  of  Theorem  4, 
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the  coefficients  a — 0>1>2,  and  a — 1)2,  are  uniquely 

determined.  If  ei,i  is  non-degenerate  at v ; , then  in  addition  is  uniquely 
determined.  Otherwise,  this  coefficient  is  not  determined.  Hence,  if  e^i  is  non- 
degenerate at  tq,  then  we  have  to  choose  exactly  one  additional  coefficient  to 
determine  all  coefficients  of  C2,  and  otherwise,  we  have  to  choose  exactly  two 
additional  coefficients.  We  claim  that  in  the  latter  case  every  choice  of  exactly 
two  additional  coefficients  from  the  set  {a£j]0,  a'Tjj0,  ag,,  fails  to 

determine  all  coefficients  of  C2. 


Proof:  Suppose  that  e^i  is  degenerate  at  tq  and  choose,  for  example,  aLT!i, 

- JTi  1 A1’U 


u3,0,0> 


«3 


a[T41 
“2,0,1  j 


(24  = 


and  For  simplicity,  we  set  ai  = al^J 0,  a2 

°i,M»  °5  = a2,i!o>  °6  = an<l  assume  that  the  remaining  coefficients  in 

C 2 are  zero.  By  Theorem  2, 


“3  ~ a4  = (1_^)a3>  «i  = (!  - |)az  + |a3, 

s $ 

a6  = (1  - 7>)a5  + aa4)  05  = (1  - -)«2,  0 = ((~r)(~  + -)  + l)ai  + -a6. 

P P 7 0 a a 

Eliminating  ajt  j <2  {3,4,5},  yields  ax  = (1  + <5^)a2,  o6  = (1  - ®)(1  + 

^fZ)a2-  By  some  elementary  computations,  we  obtain  for  the  determinant 
D of  the  corresponding  system 


D _ (~l)(r(m7  + g)  -57)2 
mrS^f2 


and  it  is  easy  to  verify  that  D — 0 iff  i is  degenerate  at  v\.  Other  choices  of 
exactly  two  additional  coefficients  from  the  set  {a^j0,  a^j|0,  a^\, 
can  be  examined  in  the  same  way,  which  proves  our  claim.  □ 

Note  that  if  e^i  is  non-degenerate  at  tq,  then  every  choice  of  exactly  one 

additional  coefficient  in  the  set  {a^j0,  a^J0,  determines  all  coefficients 

in  C2. 

We  finally  discuss  some  numerical  aspects  of  our  scheme.  A method  for 
constructing  nested  polygon  triangulations  A of  given  points  in  the  plane 
which  satisfy  the  conditions  of  Theorem  4 was  developed  in  [11],  Our  numer- 
ical tests  show  that  in  order  to  obtain  good  approximations,  it  is  necessary  to 
subdivide  some  of  the  triangles  (for  details  see  [10,11]).  Meanwhile,  we  have 
computed  such  examples  with  a high  number  of  interpolation  conditions.  We 
only  mention  here  that,  for  example,  Lagrange  respectively  Hermite  interpo- 
lation of  Franke’s  test  function  by  cubic  C^-splines  with  118,822  interpolation 
conditions  yields  an  error  of  4.66902  * 10-6  in  the  uniform  norm. 
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